Magnetism Localization in Spin-Polarized One-Dimensional Anderson— Hubbard 
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In order to study an interplay of disorder, correlation, and spin imbalance on antiferromagnetism, 
we systematically explore the ground state of one-dimensional spin-imbalanced Anderson-Hubbard 
model by using the density-matrix renormalization group method. We find that disorders local- 
ize the antiferromagnetic spin density wave induced by imbalanced fermions and the increase of 
the disorder magnitude shrinks the areas of the localized antiferromagnetized regions. Moreover, 
the antiferromagnetism finally disappears above a large disorder. These behaviors are observable 
in atomic Fermi gases loaded on optical lattices and disordered strongly-correlated chains under 
magnetic field. 

PACS numbers: 71.10.Fd, 71.10.Pm, 71.23.-k, 03.75.Ss 



An atomic Fermi gas loaded on an optical lattice 
(FGOL) has been one of the most active target in the 
atomic gas field fl*] since the successful observation of 
the superfluid-insulator transition in the Bose counter- 
part In FGOL, its interaction tunability associated 
with the Feshbach resonance and lattice formation flexi- 
bility due to optical operations offer a great opportunity 
to systematically study the Hubbard model [J| and its 
extended versions 

Thus, FGOL is expected to be 
a promising experimental reality to resolve a wide range 
of controversial issues in condensed matter physics . 

The pairing mechanism in High-Tc superconductors 
has been one of the mostly debated issues for the latest 
two decades. The superconductivity emerges by doping 
holes via chemical substitutions on the Mott-insulator 
mother phase f7|. This drastically complicates theoret- 
ical studies because the pairing occurs on the strongly- 
correlated stage influenced by disorders. This fact stimu- 
lates systematic studies for the Anderson-Hubbard (AH) 
model d, S, US IHI- FGOL is an good experimental 
testbed to solve the complicated problem [1, Q . 

In addition to such a practical motivation from the 
condensed matter side, FGOL inspires its own unique in- 
terests which lead to new frontiers. Its typical example 
is an arbitrary spin imbalance, which has been the most 
intensive target for a recent few years in atomic Fermi 
gases Therefore, in this paper, we investigate the 

spin imbalance effects on the AH model. This will be one 
of the most fruitful target as well as the balanced case 
for the future FGOL experiments. Moreover, we point 
out that its weak imbalance case has a real counterpart, 
which corresponds to disordered strongly-correlated ma- 
terials under magnetic field [isj . Thus, we examine the 
model from the weak imbalance side. 

The spin imbalance is believed to bring about a spa- 
tially inhomogeneous pairing like Fulde-Ferrel-Larkin- 
Ovchinikov state 14\ in attractive Hubbard model [3]. 



This is a widely spread idea, while the imbalance effects 
have been very little investigated in the repulsive case. 
We therefore concentrate on antiferromagnetic phases 
characteristic to the repulsive Hubbard model and study 
how not only the imbalance but also the disorder mag- 
nitude affect the phases. In this paper, we focus on 
one-dimensional (ID) repulsive polarized AH model at 
the half-filling by using the density matrix renormal- 
ization group (DMRG) method [H, [i3|. The highfight 
in this paper is that disorders localize antiferromagnetic 
phases accompanied by a phase separation from the non- 
magnetized one and finally eliminate the antiferromag- 
netism in a strong disorder range. These results are di- 
rectly observable in FGOL's by using the atomic-density 
profile probe which is presently the most standard tech- 
nique (for the experimental setup, see Ref. [5]). 
The Hamiltonian of the ID AH model is given by 
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where refers to the nearest neighbors i and j — i±l, 
t is the hopping parameter between the nearest neighbor 
lattice sites, U is the on-site repulsion, c^^ (c^J is the 
annihilation- (creation-) operator and no-i(= c^^ri'^ai) 
the site density operator with spin index a =t, |, and the 
random on-site potential is chosen by a box probability 
distribution 7'(e,) = 61(1^/2 - je^D/M^, where e{x) is the 
step function and the parameter W the disorder strength 
magnitude. In DMRG calculations, the number of states 
kept (m) is 500-700 and these numbers (m > 500) are 
enough for the most cases. In the calculations, we apply 
the open boundary condition except for the comparison 
with the periodical condition and measure the site matter 
and spin density of fermions as n-f^ ± n^i. 

Let us show DMRG calculation results on the spin im- 
balanced AH model (H]) . Figure [1] displays a random- 
ness amplitude {W/t) dependence of the on-site mat- 
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ter and spin density profile for the number of the total 
sites L = 100, the number of the spin-up and spin-down 
fermions (iV|,iV|) = (51,49) (half-filling), and U/t = 10. 
In this case, two up-spin fermions do not have their 
(down-spin) partners. The matter density profile is al- 
most flat for W < U a,s seen in Figs. [IIa)-[lJc), while 
it is drastically disturbed for W > U. This is com- 
pletely the same as that of the well-known balanced case. 
However, the spin density profile is significantly differ- 
ent. Firstly, in the clean case as shown in Fig. [TJa) 
{W/t = 0), one finds an antiferromagnetic spin density 
wave (ASDW) whose periodicity is found to be inversely 
proportional to — Ni [e.g., see Figs. [3ljh) and ^i), 
in which — Ni = 2 and 4, respectively]. Here, it is 
noted that any ASDW phases are never observed in the 
perfectly balanced cases irrespective of the presence of 
randomness [l^l- This clearly indicates that the imbal- 
ance is responsible for the ASDW phase. Secondly, as 
one increases the disorder strength, a part of the ASDW 
(amplitude) "locally" vanishes and the depressed regions 
expand as seen in Figs, [ijb) and[ljc). This tendency 
becomes remarkable when W exceeds over U as seen in 
Figs. [Hd) and[T]^e), in which two ASDW phases are lo- 
calized and isolated each other. This isolation can be 
explained by the complete localization of two excess up- 
spin fermions for W > U since the ASDW can be created 
only at the localized spots of the excess up-spin particles. 
In fact, the change in the up and down spin profiles in 
Fig- [2c)"{IIe) supports the idea. The further increase of 
W diminishes the antiferromagnetic structure as seen in 
Fig.[ljf), in which the localized structure is characterized 
by positive peaks instead of the staggered (plus-minus) 
moment alternation. This is because the strong random- 
ness fully dominates over the other effects. However, we 
note here that it is difficult in this large disorder range 
to judge whether the result [Fig.^f)] is the true ground- 
state or not. The reason is that in this strong W range 
tiny changes of W (e.g., W + dW) give entirely differ- 
ent spin-density distributions in non-continuous manner, 
which is not observed when < W/t < 14. Generally, 
it is well-known in strongly glassy situations that a tiny 
change in calculation parameters results in a drastic dif- 
ferent consequence. Thus, we expect that there are a lot 
of local minima in this strong disorder range. Fig. [Ijf ) is 
a localization profile selected among those minima. 

In order to characterize the spin density profile in a 
wide range of U and W, we define the following function. 



S{U,W) = {J2 



\n^,{e,U,W)-ni,{e,U,W)\ 
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where nai{e, U, W) is the local site density under a ran- 
dom potential symbolically specified by e at a certain set 
of U and W, and means an algebraic average for 
various random realizations. From the expression, it is 
found that S{U, W) gives an indicator how large the spin 



(a)TI'/t = T + i^ T-i 





A'AV>;»■v;J•./••^^\^A»Vtf^'>^■^"r■,^y..«Al/ltt 

20 40 60 80 100 20 40 60 80 100 



FIG. 1: The randomness magnitude W dependence of the 
matter and spin density profiles (n|i±n|i respectively) at the 
half-filling for U/t = 10 with (A^T.A^'i) = (51,49). A profile 
of the selected random potential is depicted on the bottom of 
each figure in an arbitrary unit (grey dashed line). In figures 
(c)-(e), up and down spin density profiles are shown. For all 
calculations, m is 500. 



moment develops on each site. If the staggered moment 
widely grows, then S{U, W) gives a relatively large value. 
Thus, a map of S{U, W) in a wide range of U and W is ex- 
pected to clarify an interplay of U and W on the ASDW 
phase localization. Figure [UJa) shows a contour plot of 
S{U,W), which is averaged over ten realizations of ran- 
dom potentials for L = 100 and (A^t-^i) = (51,49). In 
this figure, when one increases W along a fixed U/t line 
(e.g., U/t = 10 hue) from W/t = to 20, it is found that 
the averaged moment of ASDW very slowly decreases in- 
side the region W < U. This is consistent with the spin 
density profile as seen in Figs. ma)-[TJc) in which the 
areas of ASDW phases slowly diminish with increasing 
W. When W exceeds over U, the variation of S{U, W) 
suddenly changes to a fast suppression. This reflects the 
change in the localization of the excess up-spin particles 
as seen in Figs.[Hc) and[Ud). 

In order to qualify the present randomness averaging 
on Fig.[2Ka), we introduce the following function 



D{U, W) 



{[s,{u,w)/s{u,w)]-iy 



(3) 



This corresponds to a standard deviation on the random- 
ness average. Fig. ^h) is a contour map of D{U, W) in 
the same range as Fig.[2l^a). One finds that in a small W/t 
range D{U, W) shows a very small value (almost zero) 
and D{U, W) reaches about 0.2 around W = U + 3t At 
line when increasing W/t at a constant U/t. Thus, the 
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FIG. 2: Contour plots of the value of (a) S{U, W) [Eq. (O] 
and (b) D{U,W) [Eq. (|3}] with L = 100 and (A^T.^i) = 
(101, 99) at the half-fiUing in a range oiU {Q <U < 20) and 
W {0 <W <2Q). The step value for both U/t and W/t is 2. 
For all calculations, m is 500. 



averaged values in Fig. [Ha) are sufficiently qualified ex- 
cept for above W + it line. This result means that 
the qualitative change as observed around W = U va 
Fig. [2ja) is not a side effect associated with the averag- 
ing but an essential feature in this system. 

Next, we investigate the polarization strength depen- 
dence of the spin density profile. Figures [3Ua)-[3tg) dis- 
play the correspondent results made at the half-filling in 
L = 200, U/t = 18 and W/t = 20. In the range of 
these parameters, clear localization of the ASDW phases 
can be observed with the phase separation from the non- 
magnetized phases in a slight polarized case [e.g., see 
Fig. [Hd)]. Firstly, Fig. [3l^a) displays the charge and 
spin density distributions in {N-^,N^) — (101,99). One 
finds two magnetized regions, in ASDW phases are lo- 
calized with the localization of two extra fermions. The 
slight increase of the polarization [(A'^i, A^j) (102,98)] 
increases the number of the magnetized regions as seen 
in Fig. [31(b). Here, we note that the magnetized regions 
formed in the less polarized case as Fig. [3ja) are kept. 
This is in contrast to the clean systems {W/t = 0) as 
shown in Figs, ^h) and^i) whose polarizations are the 
same as Figs.[3lja) and[3ljb), respectively. This is a typi- 
cal feature characteristic to disordered systems, in which 
memory effects can be frequently observed. The further 
increase of the polarization extends the magnetized re- 
gions as shown in Figs. [3Kc)-[3Kg). In Fig.[31[g), the mag- 
netized (ASDW) regions cover all sites, which are almost 
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FIG. 3: The polarization dependence of the spin density pro- 
files [(a)-(g)] with U/t = W/t = 20 and that without random 
potential [(h) and (i)] at U/t = 20. The last two figures (h) 
and (i) without the randomness are displayed for a compari- 
son with disordered cases (a) and (b) respectively. The num- 
ber of up and down spin particles {Ni,Ni) are (a) (101,99), 
(b) (102,98), (c) (103,97), (d) (104,96), (e) (105,95), (f) 
(110,90), (g) (120,80), (h) (101,99), and (i) (102,98), re- 
spectively. The random potential shape are displayed on the 
bottom of the figures (a)-(g) (gray dashed hues). For all cal- 
culations, m is 600. 

positively polarized although the alternation profile still 
remains. 

We investigate a system size dependence of this mag- 
netism localization to confirm that it is not a small size 
effect. In Figures [3Ka)-[Hc), we show the charge and spin 
density distributions at U/t = W/t = 20 in i = 100, 
200, and 300 with (N^.N^) = (51,49), (102,98), and 
(153,147), respectively. We clearly find that these all 
cases exhibit the same qualitative behavior, i.e., the mag- 
netized regions are localized with the separation from the 
non-magnetized regions. From these figures, we find that 
the observed magnetism localization is an intrinsic effect. 

Finally, we re-examine the model (|T]) under the peri- 
odic boundary condition to check the effect of the bound- 
ary condition. Figures [ni^a)-[51Jf) show the W dependent 
charge and spin density profiles in the periodic condi- 
tion. In Fig. [5l^a) {W/t = 0), we find complete fiat dis- 
tributions in both the charge and spin densities. They 
are characteristic to the periodic boundary condition in 
which the excess fermions fully distribute homogeneously. 
This is in contrast to the open boundary condition [com- 
pare it with Fig. [ija)]. When the randomness is added 
into the system, the ASDW phases are induced as seen 
in Figs.[n);b) undM^c). Thus, one finds that the ASDW 
phase requires two conditions, i.e., the imbalance and the 
translational symmetry breaking. When the randomness 
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FIG. 4: The size dependence of the matter and spin density 
profiles. The same polarization ratio is kept, (a) L = fOO case 
with (iVt.^i) = (51,49), (b) L = 200 case with {N^.N^) = 
(102,98), and (c) L = 300 case with (A^t.^i) = (153,147), 
respectively. U/t = W/t — 20 in all cases. The random 
potential profiles are shown on the bottom of each figure (grey 
dashed lines). For all calculations, m is 600. 
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FIG. 5; The randomness W dependent matter and spin den- 
sity profiles (n|i±n|i respectively) with L — 100, {Ni,Ni) = 
(51,49), and U/t — 10, under a random potential depicted 
on the bottom of each figure in arbitrary unit (grey dashed 
line). The periodic boundary condition is employed. For all 
calculations, m is 700. 



strength increases, one finds that the amplitude of the 
ASDW increases [see Figs. [Sfb) andj^c)]. This implies 
that the localization of extra two fermions proceeds with 
increasing W/t. At W/t = 12 and 14 [see Figs.OJd) and 
[Sje) for W > U], we find the phase separation from the 
non-magnetized phases, which is the same as the open 
boundary case. Further increase of W/t brings about 
more tight localization [Fig. OJe)] and disappearance of 
the staggered moment profile [Fig. [5{f)]. This is also the 
same as the open boundary case. 

In conclusion, we systematically studied the polarized 
AH model at the half-filling and found that the disor- 
der localizes the ASDW phases induced by the excess 
fermions. As the randomness strength increases, the ar- 
eas of the localized ASDW phases shrink with the ex- 
pansion of the non-magnetized areas, and the antiferro- 
magnetism finally vanishes. These novel disorder effects 
on polarized strongly-correlated systems are observable 
in not only ID FGOL's Q but also strongly-correlated 
disordered chains under the magnetic field. 
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